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Abstract 

We extend a recent calculation of the nuclear energy density functional in the sys- 
' tematic framework of chiral perturbation theory by computing the isovector spin-orbit 

terms: (V pp — V pn) • {Jp — Xi) Gso{kf) + {Jp — Jn)'^ Gj{kf)- The calculation includes the 
I one-pion exchange Fock diagram and the iterated one-pion exchange Hartree and Fock 

Q ' diagrams. From these few leading order contributions in the small momentum expan- 

sion one obtains already a good equation of state of isospin-symmetric nuclear matter. 
Q I We find that the parameterfree results for the (density-dependent) strength functions 

' Gsoikf) and Gj{kf) agree fairly well with that of phenomenological Skyrme forces for 

densities p > po/W. At very low densities a strong variation of the strength functions 
Gso{kf) and Gj{kf) with density sets in. This has to do with chiral singularities 'rn~^ 
\ and the presence of two competing small mass scales kf and m^. The novel density 

dependencies of Gso{kf) and Gj{kf) as predicted by our parameterfree (leading order) 
. calculation should be examined in nuclear structure calculations. 

(N 



o 



PACS: 12.38.Bx, 21.30.Fe, 21.60.-n, 31.15.Ew 



^ ' Among the various phenomenological interactions that have been used extensively in the 

C ■ description of nuclei, the Skyrme force [1] has gained much popularity because of its analytical 
simplicity and its ability to reproduce nuclear properties over the whole periodic table within 
^ i the self-consistent Hartree-Fock approximation. Several Skyrme parameterizations have been 
I tailored to account for single-particle spectra [2], giant monopole resonances [3] or fission 
barriers of heavy nuclei [4]. Recently, a new Skyrme force which also reproduces the equation 
of state of pure neutron matter up to neutron star densities, pn — 1.5 fm~^, has been proposed 
in ref. [5] for the study of nuclei far from stability. A microscopic interpretation of the various 
parameters entering the effective Skyrme forces is generally put aside. Sometimes the energy 
density functional is just parameterized without reference to any effective (zero-range) NN- 
interaction [6]. 

Another widely and successfully used approach to nuclear structure calculations are rela- 
tivistic mean- field models [7, 8] . In these models the nucleus is described as a collection of in- 
dependent Dirac-particles moving in self-consistently generated scalar and vector mean-fields. 
The footprints of relativity become visible through the large nuclear spin-orbit interaction 
which emerges in that framework naturally from the interplay of the two strong and counter- 
acting (scalar and vector) mean-fields. The corresponding many-body calculations are usually 
carried out in the Hartree approximation, ignoring the negative-energy Dirac-sea. For a recent 
review on self-consistent mean-field models for nuclear structure, see ref. [6]. In that article 
the relationship between the relativistic mean-field models and the Skyrme phenomenology is 
also discussed. 
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The first necessary conditions to be fulfilled by any phenomenological nuclcon-nucleon 
interaction come from the (few empirically known) properties of infinite nuclear matter. These 
are the saturation density po = 2Ar^o/37r^5 the binding energy per particle —E{kfo) and the 
compression modulus K — /cjo-E"(A;/o) of isospin-symmetric nuclear matter as well as the 
asymmetry energy A{kfo). In a recent work [9], we have applied the systematic framework of 
chiral perturbation theory to the nuclear matter many-body problem. In this calculation the 
contributions to the energy per particle, E{kf), originate exclusively from one- and two-pion 
exchange between nucleons and they are ordered in powers of the Fermi- momentum kf (modulo 
functions of kf/m.„). It has been demonstrated in ref.[9] that the empirical saturation point 
(po — 0.17 fm~^ , E{kfo) ~ —15.3 MeV) and the nuclear matter compressibility K ~ 255 MeV 
can be well reproduced at order 0{kj) in the small momentum expansion with just one single 
momentum cut-off scale of A ~ 0.65 GeV which parameterizes the short-range NN-dynamics 
necessary for nuclear binding. Most surprisingly, the prediction for the asymmetry energy, 
A{kfo) = 33.8 MeV [9], is in good agreement with its empirical value. In fact very similar 
results for nuclear matter can be obtained already at order 0{kj) in the small momentum 
expansion (by dropping the rclativistic l/A/^-corrcction to Ivr-exchange and the irreducible 
27r-exchange of order 0{k^)) with a somewhat reduced cut-off scale of A ~ 0.61 GeV (for 
detailed results see ref.[10]). 

Given the fact that many properties of nuclear matter can be well described by chiral irN- 
dynamics treated perturbatively up to three- loop order it is natural to consider in a further step 
the energy density functional relevant for inhomogeneous many-nucleon systems (i.e. finite 
nuclei). Such an extension to inhomogeneous many-nucleon systems can be done with the 
help of the density matrix- expansion of Negele and Vautherin [11]. The bilocal density- matrix 
(given by a sum over the occupied energy eigenfunctions) is expanded in relative and center- 
of-mass coordinates with expansion coefficients determined by purely local quantities: nucleon 
density p(r), kinetic energy density T(r ) and spin-orbit density J(r ). The Fourier-transform 
of the (so expanded) density-matrix defines in momentum-space a "medium-insertion" T{p,q) 
for the inhomogeneous many-nucleon system which then allows to calculate diagrammatically 
the nuclear energy density functional. 

In a recent work [12] we have considered the isospin-symmetric case of equal (even) proton 
and neutron number. The corresponding energy density functional takes the general form: 
S[p, r, J] = pEikf) + [r - 3pk}/5]/2M*{p) + {Vpf F^ikj) + Vp ■ JFUkf) + P Fj{kf). We 
have found that the effective nucleon mass M*{p) deviates at most by ±15% from its free space 
value M, with 0.89M < M*{p) < M for p < 0.11 fm"^ and M*{p) > M for higher densities 
P < Po — 0.174 fm~^. Interestingly, a recent large scale fit of (almost two thousand^ nuchde 
masses by Pearson et al. [13] finds a similarly enhanced effective nucleon mass: M*(po) ~ 
1.05M. The strength of the (Vp)^-term, Fv(A;/o), is (at saturation density) comparable to 
that of phenomenological Skyrme forces. The magnitude of Fj{kfo) accompanying the squared 
spin-orbit density J ^ comes out considerably larger. Both quantities increase strongly as the 
nucleon density p tends to zero. This has to do with the explicit presence of the small mass 
scale m^r = 135 MeV which amounts to about half of the Fermi momentum kfo in equilibrated 
nuclear matter. The strength of the isoscalar nuclear spin-orbit interaction, Fso{kfo), as given 
by iterated Ivr-exchange is (at saturation density) about half as large as the corresponding 
empirical value ~ 90 MeVfm^, however, with the wrong negative sign. This isoscalar spin-orbit 
interaction is not a rclativistic effect but proportional to the large nucleon M. From that result 
it becomes clear that perturbative chiral pion-nucleon dynamics (alone) cannot account for 
the mechanisms underlying the empirical isoscalar nuclear spin-orbit force, whereas relativistic 
scalar- vector mean-field models [7, 8] give a successful phenomenology of it. 
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The purpose of the present work is to calculate using the same framework as in ref. [12] 
the isovector spin-orbit interactions generated by one- and two-pion exchange in order to 
reveal the underlying isospin dependence. Note, for example, that the Skyrme force gives 
rise to isoscalar and isovector spin-orbit interactions with a fixed ratio 3:1. In particular, this 
restrictive feature of the Skyrme energy density functional has been made responsible for the 
less accurate description of isotope shifts in the Pb region [14] in comparison to relativistic 
mean-field calculations. Naturally, one expects that the finite-range character of the two-pion 
exchange interaction will lift such restrictive features of the zero-range Skyrme force. 

Let us begin with writing down the explicit form of the isovector spin-orbit terms in the 
nuclear energy density functional: 

£[pp, Pn, Jp, Jn] = (Vpp - Vpn) ' {Jp - Jn) Gso{kf) + {Jp - Jnf Gj{kf) + ... (l) 

Here, 

Pp^^) - %P - E '^lfir)¥;:lir) , (2) 

aeocc 

are the local proton and neutron densities written in terms of the corresponding (local) proton 
and neutron Fermi- momenta A;p „(r ), and expressed as sums over the occupied single-particle 
orbitals ^g(f ). The spin-orbit densities of protons and neutrons are defined similarly: 

JpAr) = E Hfim ^ X V*g(f ) . (3) 

a&occ 

Furthermore, Gso{kf) and Gj{kf) denote the associated strength functions. In Skyrme param- 
eterizations [2, 3, 4, 5, 13] these are just constants, Gso{kf) = Wq/A and Gj{kf) = {ti — 12)/32, 
whereas in our calculation their explicit density dependence originates from the finite-range 
character of the Itt- and 27r-exchange interaction. 

The starting point for the construction of an explicit nuclear energy density functional 
S[pp, Pn, Jpj Jn] is the bilocal density-matrix as given by a sum over the occupied energy eigen- 
functions: J2aeocc^^p'l{''^ ~ ^/2)^p'ri^(^+ a/2). According to Negele and Vautherin [11] it can 
be expanded in relative and center-of-mass coordinates, a and f, with expansion coefficients 
determined by local quantities (nucleon density, kinetic energy density, spin-orbit density). 
As outlined in section 2 of ref. [12] the Fourier-transform of the (so expanded) density- matrix 
defines in momentum-space a medium-insertion r{p,q) for the inhomogeneous many-nucleon 
system. It is straightforward to generalize this construction to the isospin-asymmetric sit- 
uation of different proton and neutron local densities Pp,n{^) and Jp^f ). We display here 
only that part of the medium-insertion r(p, q) which is actually relevant for the diagrammatic 
calculation of the isovector spin-orbit terms eq.(l): 

r(p,g) = j dhe-^^-^^\^-^e{kp-\p\) + ^-^e{kr.-\p\) 

na-\px{Jp- Jn)\ + ...y (4) 

When working to quadratic order in deviations from isospin symmetry (proton-neutron differ- 
ences) it is sufficient to use an average Fermi- momentum kf in. the prefactor of the isovector 
spin-orbit density Jp — Jn. The double line in the left picture of Fig. 1 symbolizes this medium 
insertion together with the assignment of the out- and in-going nucleon momenta p±g'/2. The 
momentum transfer q is provided by the Fourier components of the inhomogeneous (matter) 
distributions Pp,n{^) and Jp^f ). 
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Fig. 1: Left: The double line symbolizes the medium insertion defined by eq.(4)- Next are 
shown: The two-loop one-pion exchange Fock-diagram and the three-loop iterated one-pion 
exchange Hartree- and Fock- diagrams. Their combinatoric factors are 1/2, 1/4 and 1/4, in 
the order shown. 

Now we turn to the analytical evaluation of the pion-exchange diagrams shown in Fig. 1. 
We give for each diagram only the final result omitting all technical details related to extensive 
algebraic manipulations and solving elementary integrals. A collection of the relevant "master 
integrals" can be found in the appendix of ref.[12]. We obtain from the l7r-exchange Fock 
diagram with two medium insertions: 



Gj{kf) 



9\ 



10 + 24^2 



(8m,/0n 3(1 + 4x^2)^ 



u 



(5) 



This contribution to Gj{kf) is just —1/3 of the contribution to the isoscalar strength func- 
tion Fj{kf) (see eq.(lO) in ref. [12]) as a simple consequence of the isospin-trace tr[Ta{Js + 
'TzJv)'Ta{Js + Ta^-u)] = 6 — 2 J^. There is no contribution of the Ivr-exchange Fock-diagram to 
the strength function Gso{kf) since the spin-trace with one insertion proportional to <?■ (px J„) 
vanishes. We note also that the l7r-exchange Hartree diagram (not shown in Fig. 1) van- 
ishes identically either because of a zero spin-trace or because of the momentum integral 
/ d^ppd{kp^n ~ = 0. The iterated one-pion exchange Hartree diagram with two medium 
insertions leads to the result: 



(6) 



which is —2/3 of the contribution to the isoscalar spin-orbit strength Fgoikf) (see eq.(13) in 
ref.[12]). Let us briefly explain the mechanism which generates the strength function Gso{kf). 
The exchanged pion-pair transfers the momentum g*between the left and tlie right nucleon ring 
and this momentum q enters also the pseudovector vriV-interaction vertices. The isovector spin- 
orbit strength Gso{kf) arises from the spin-trace tr[(f-{l+q/2) (T-{l—q/2) (T-{px J^)] = 2i{qxl)- 
{px Jy) where q gets converted to V/Cp — VA;„ ~ 7r^(Vpp — Vp^) / kj by Fourier transformation. 
From the iterated one-pion exchange Fock diagram with two medium insertions we obtain: 



Gsoikf) 



5g\M 



67rm, (4/^)4 
1 . 14 



arctan 2u 



3 + 



In 



u 



2 



arctan u 



2u^ 1 + 4^2 



2m2) u{1 + 2u^j 

arctan 2x — arctan x 
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7/2(1 + 2x2) 
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Gj{kf) 



bg\M 



37rm^ (8/^)4 U(l + 2^2) 



2(5 + 8^2) 2 arctan 2m 

arctan u _ + 



1 1 + 



+^ln 



dx 



arctan x — arctan 2x 
1 + 2x2 



(8) 



One notices a relative factor of —5/3 in comparison to the contributions to the isoscalar 
strength functions Fso,j{kf) (see eqs.(16,17) in ref.[12]) which comes from the isospin-trace 
ii[Ta{Js + T'iJv)TbTa{Js + T'iJv)Tb] — lOJ^ — 6J5 of that diagram. In case of the iterated 
one-pion exchange Hartree diagram with three medium insertions one has to evaluate nine- 
dimensional principal value integrals over the product of three Fermi spheres. We find the 
following contribution to the isovector spin-orbit strength: 



37r2m^(4/^) 



dxx 



1 j' 2Mg4[2xV(g - s') + u^{2s' - s)] 

-1 I (1 + S2)2(t,2 _ ^2^2) 



+ 



u{3u - hx^y ) 



^2^2 



AxyH 



3 arctan s 



3s + 2s^ 



l + s2 



{8xy - 5s - s^)s'^ - 2xys'{7s + 3s^) 



(1 + S2)3 

+s'^xy{ll + 7s^) + (s + s^){2xy - s)(s" - s') 



a;2|/2. 



(9) 
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with the auxiliary functions H = \n(u + xy) — hi{u — xy) and s = xy -\- \Ju^ — x'^ 
partial derivatives of s are abbreviated by s' = uds/du and s" — u^d^s/du^. In comparison 
to the contribution to the isoscalar spin-orbit strength Fso{kf) (see eq.(20) in ref.[12]) relative 
isospin factors of ±2/3 have occurred in various subparts. The contribution of the iterated 
one-pion exchange Hartree diagram with three medium insertions to the strength function 
Gj{kf) reads: 



Gjikf) 



9aM 



37r2m^(4A)4 



1 + 12^2 - 24?/ - 96^6 1 + 2^2 



iu^y'^ 



Jo ^^(1+ 4^2^2)4 



M(l + 4ti2)3 41*3 

5 - (30 + ?,2u'^)y'^ 

l + y_ 



ln(l + 4^2) 



+ (35 + 24^2 - 16M4)y4 + 56k^/ + 48^^^^ 



In 



(10) 



It is obtained if only if both insertions proportional to a ■ {'pi,2 x ) stand under a single 
spin-trace and this feature explains also the relative isospin factor —1/3 in comparison to the 
contribution to Fj(kf) written in eq.(21) of ref.[12]. The iterated one-pion exchange Fock 
diagram with three medium insertions is most tedious to evaluate. It is advisable to split the 
contributions to the strength functions Gso,j{kf) into "factorizable" and "non-factorizable" 
parts. These two pieces are distinguished by whether the nucleon propagator in the denom- 
inator can be canceled or not by terms from the product of TriV-interaction vertices in the 
numerator. We find the following "factorizable" contributions: 
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Gj{kf) = 
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(12) 



+4Q£dx I [2 + 2m^ - 3(1 + 'u^)^^-^ - 3x^]uL'^ + 2[3x-^{u^ + u"^) - 2 - 



+ (1 + UX- u^)[l + {u + a;)X + {I - ux - u^)[l + {u - xy] 



21-1 



L - 3u^x 



3^-2 



with the auxihary function: 



1 , 1 + iu + xY 
L — — m 



Ax l + {u-xy ^^^^ 

The non-factorizablc contributions from the iterated one-pion exchange Fock diagram with 
three medium insertions (stemming from nine-dimensional principal value integrals over the 
product of three Fermi spheres) read on the other hand: 



Gsoikf) 
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(14) 



Gj{kf) 
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Au'z' - ln(l + Au'z') 
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m(1 + 4m2i/2)3 
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-I- / dr 

lo 4^10(1 + S2)2(l+t2)2 



(1 + 4m2|/2)2(1 + 4^^2^2) 

[s + s'){s"-s') + {3 + s')s''^ 



X 



(i + i3)(i"-i') + (3 + i')i" 



(15) 



with the auxiliary function t = xz + \Jv? — x'^ ^ x^z'^ and its partial derivatives i' = udt/du 
and t" — u^dH/du^. In comparison to the contributions to the isoscalar strength functions 
Fso,j{kf) (see eqs. (25,26,29,30) in ref.[12]) there occur relative isospin factors of —1/3 and —5/3 
in various subparts. Let us add some general power counting considerations for the nuclear 
energy density functional S[pp, pn, Jp, Jn]- Counting the Fermi momenta kp^nj, the pion mass 
m-jr and a spatial gradient V collectively as small momenta one deduces from eqs. (2,3) that 
the nucleon densities Pp,n(^) and the spin-orbit densities Jp,n(^) are quantities of third and 
fourth order in small momenta, respectively. With these counting rules the contributions from 
iTT-exchange to the nuclear energy density functional S[pp, Pn, Jp, Jn] are of sixth order in small 
momenta while all contributions from iterated Ivr-exchange are of seventh order. Concerning 
NN-interactions induced by pion-exchange the isovector spin-orbit terms presented here are 
in fact complete up-to-and-including seventh order in small momenta. 
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Next we turn to the numerical results. We use the physical parameters: M — 939 McV 
(nucleon mass), m^r = 135 MeV (neutral pion mass), /^r = 92.4 MeV (pion decay constant) 
and gA = 1-3 (equivalent to a vriV-coupling constant Qt^n = QaM/ f^^ = 13.2). The full line 
in Fig. 2 shows the result of iterated l7r-exchange for the strength function Gso{kf) belonging 
to the isovector spin-orbit coupling term (Vpp — Vpn) • {Jp — Jn) versus the nucleon density 
p = 2kj/3n'^. For comparison we have drawn the constant values Gso{kf) — Wq/A of the 
three Skyrme forces Sly [5], MSk [13] and SkP [15] (horizontal dashed lines). In the case of 
Sly and MSk we have performed averages over the various parameter sets Sly4-7 and MSkl-6. 
One observes a fair agreement of our parameterfree prediction with these empirical values. In 
contrast to the isoscalar spin-orbit strength Fso{kf) (see Fig. 4 in ref.[12]) iterated l7r-exchange 
gives now the correct positive sign. The density dependence of the isovector spin-orbit strength 
Gsoikf) is moderate for densities p > po/10). A rapid decrease sets however in when p tends 
to zero. 
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0.00 0.05 0.10 0.15 0.20 

P [fm '] 

Fig. 2: The strength function Gso{kf ) related to the isovector spin-orbit term {Vpp — Vpn) ■ 
{Jp — Jn) in the nuclear energy density functional versus the nucleon density p = 2k'f/37i^. The 
three horizontal dashed lines show the constant values Gsoikf) = Wq/A of the Skyrme forces 
Sly [5], MSk [13] and SkP [15]. 

Finally, we show in Fig. 3 the strength function Gj{kf) accompanying the squared isovector 
spin-orbit density [Jp — Jnf' in the nuclear energy density functional versus the nucleon density 
p — 2A;y/37r^. For comparison we have drawn the constant values Gj{kf) = [ti — t2)/32 of 
the three Skyrme forces Sly [5], SkP [5] and MSk [13], (dashed lines). One observes that our 
parameterfree prediction for the strength function Gj{kf) lies well within the band spanned 
by these empirical values. Again, the density dependence of Gj{kf ) is moderate for densities 
p > po/5 and a rapid decrease sets in when p tends to zero. Although not visible each (full) 
curve in Figs. 2,3, approaches a finite (negative) value at p = 0. One can analytically derive 
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the following low density limits: 



GsoiO) 
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37rm,(4/,)4 



-33.8MeVfni^ , 



9\ 



1 + 



= -108.0 MeVfm^ 



(16) 



(17) 



2567r/2 

to which only the diagrams with two medium insertions contribute. The large numbers in 
eqs.(16,17) arise from negative powers of the pion mass (so-called chiral singularities). 
The most singular m~^-term can be traced back to the l7r-exchange Fock diagram. Note also 
the relation between the isovector and isoscalar spin-orbit strengths at zero density, Gso(O) = 
-^so(0)/3 (for Fso(O) see eq.(42) in ref.[12]). This is a necessary consistency check on our 
diagrammatic calculation. At extremely low densities [kj << 171^^/2) even the pion-cxchange 
interaction becomes effectively short-ranged and therefore the constraint Gso{kf ) = Fso{kf)/3 
known from the (zero-range) Skyrme spin-orbit force must hold. 

In this context is important to keep in mind that if pionic degrees of freedom are treated 
explicitly in the nuclear matter problem the low density limit is realized only at extremely low 
densities kf << Often, the opposite limit where the pion mass nij^ can be neglected 

against the Fermi momentum kf is already applicable at the moderate densities relevant for 
conventional nuclear physics. This is exemplified here by the approximate density dependence 
Gso,j{kf) ~ kj^ (see Figs. 2,3). Such a p~^/^-behavior becomes exact in the chiral limit m^r = 
as can be deduced by simple mass dimension counting of the iterated Ivr-exchange diagrams 
(the basic argument is that M/ f^kf has the correct unit of MeVfm^). 
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Fig. 3: The strength function Gj{kf) accompanying the squared isovector spin-orbit density 
{Jp — in the nuclear energy density functional versus the nucleon density p = 2/^^/37?^. 
The three horizontal dashed lines show the constant values Gj{kf) = (ti— 12)/32 of the Skyrme 
forces Sly [5], SkP [15] and MSk [13]. 
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In summary we extended in this work our recent calculation of the nuclear energy density 
functional [12] in the framework of chiral perturbation by computing the isovector spin-orbit 
terms: (Vpp — Vp.„) ■ {Jp — Jn) Gso(kf) + ( Jp — J„)^ Gjikf). Our calculation includes the Itt- 
exchange Fock diagram and the iterated Ivr-exchange Hartree and Fock diagrams. These few 
leading order contributions in the small momentum expansion give already a good equation of 
state of isospin symmetric infinite nuclear matter [9, 10]. The step to inhomogeneous many- 
nucleon systems is done with the help of the density-matrix expansion of Negele and Vautherin 
[11]. The specific isospin-structures of the Ivr- and 27r-exchange interaction show up through 
relative isospin-factors —1/3, ±2/3 and —5/3 of various subcontributions. We find that the 
parameterfree results for the (density-dependent) strength functions Gso{kf) and Gj{kf) agree 
fairly well with that of phenomenological Skyrme forces for densities p > po/lO. At very low 
densities a strong variation of the strength functions Gso{kf) and Gj{kf) with density sets in. 
This has to do with chiral singularities and the presence of two competing small mass 
scales kf and myr. The novel density dependencies of Gsoikf) and Gj{kf) as predicted by our 
parameterfree (leading order) calculation should be examined in nuclear structure calculations. 
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